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ABSTRACT 

One of the proposed damping mechanisms of coronal (transverse) loop oscillations in the 
kink-mode is resonant absorption as a result of the Alfven speed variation at the outer boundary 
of coronal loops. Analytical expressions for the period and damping time exist for loop models 
with thin non-uniform boundaries. They predict a linear dependency of the ratio of the damping 
time to the period on the thickness of the non-uniform boundary layer. Ruderman and Roberts 
used a sinusoidal variation of the density in the non-uniform boundary layer and obtained the 
corresponding analytical expression for the damping time. Here we measure the thickness of 
the non- uniform layer in oscillating loops for 11 events, by forward- fitting of the cross-sectional 
density profile n e (r) and line-of-sight integration to the cross-sectional fluxes F(r) observed 
with TRACE 171 A. This way we model the internal rij and external electron density n e of the 
coronal plasma in oscillating loops. This allows us to test the theoretically predicted damping 
rates for thin boundaries as function of the density ratio \ — n e /rii. Since the observations 
show that the loops have non-uniform density profiles we also use numerical results for damping 
rates to determine the value of x f° r the loops. We find that the density ratio predicted by 
the damping time, xleda = 0.53 ± 0.12, is a factor of w 1.2 — 3.5 higher than the density 
ratio estimated from the background fluxes, x — 0.30 ± 0.16. The lower densities modeled 
from the background fluxes are likely to be a consequence of the neglected hotter plasma that 
is not detected with the TRACE 171 A filter. Taking these correction into account, resonant 
absorption predicts damping times of kink-mode oscillations that are commensurable with the 
observed ones and provides a new diagnostic of the density contrast of oscillating loops. 



Subject headings: Sun : corona — Sun : magnetic fields — waves 



1. INTRODUCTION 

Oscillations of coronal loops have now been detected virtually in all wavelengths (for a recent review 
see, e.g., Aschwanden 2003). Most of these oscillations have been interpreted in terms of standing 
(eigen-modes) and propagating MHD waves (for a recent theoretical review see, e.g., Roberts & Nakariakov 
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2003). The MHD eigen- modes include fast sausage and kink modes that produce transverse oscillations 
with Alfvenic speed, slow magneto-acoustic modes that produce longitudinal oscillations with sound speed, 
as well as torsional modes that produce sheared azimuthal oscillations. Obviously, observations of such 
oscillating systems provide direct measurements of Alfven speeds and sound speeds, which in combination 
with electron density measurements can be used to infer the magnetic field in the corona, which is very 
difficult to determine by other means. This important new diagnostic has been dubbed coronal seismology 
(Roberts, Edwin, & Bcnz 1984; Roberts & Nakariakov 2003). 

Most of the coronal loops that exhibit oscillations have been found to be strongly damped, typically 
having an exponential damping time <d of a few oscillation periods P (Nakariakov et al. 1999; Schrijver 
et al. 2002; Aschwanden et al. 2002). Theoretical models of damping mechanisms include: (1) non-ideal 
effects such as viscous and Ohmic damping, optically-thin radiation, thermal conduction, (2) wave leakage 
across the sides of the loop boundaries, (3) wave leakage at the chromospheric footpoints, (4) phase mixing 
in inhomogeneous loop regions, and (5) resonant damping at the sides of loop boundaries. The first 
three effects are believed to be weak for fast kink-mode oscillations, while the latter two are considered 
as most important (Goossens 1991; Poedts 2002; Ruderman & Roberts 2002; Ofman & Aschwanden 
2002; Goossens, Andries, & Aschwanden 2002; Erdelyi 2003). First observational tests with TRACE data 
revealed that the scaling law of the damping time as function of other physical parameters (loop length 
L and period P) favors the phase mixing mechanism (Ofman & Aschwanden 2002), but the mechanism 
of resonant absorption can explain the observed damping times equally well if the inhomogencity length 
scale is a fraction of w 15 — 50% of the loop radius (Goossens et al. 2002). More accurate tests to decide 
between these two damping mechanisms require the knowledge of the inhomogeneity length scale I and the 
density ratio n e /rii between the external and internal electron density of oscillating loops. The knowledge 
of the density ratio n e /rii is also required to calculate a coronal magnetic field strength B from a loop 
with oscillation period P and length L, which is a fundamental tool of coronal seismology (Nakariakov 
& Ofman 2001). In this paper here we measure for the first time these additionally required parameters 
in 11 kink- mode oscillations events, for which the damping times have been reliably determined earlier 
(Aschwanden et al. 2002). This allows for a more rigorous quantitative test of the damping mechanisms, 
with no free parameters for the theoretical model of resonant absorption. We find that the mechanism of 
resonant absorption is commensurable with the observed damping times. The data analysis and discussion 
of observational parameters are discussed in Section 2, while conclusions are summarized in Section 3. 

2. DATA ANALYSIS 

We analyze 11 loop oscillation events from the study of Aschwanden et al. (2002) for which a reliable 
damping time to has been determined (e.g., see event #1 in Fig. 1). The same data set of these 11 events 
is also studied in Ofman & Aschwanden (2002) and in Goossens, Andries, & Aschwanden (2002). The date 
and times of the observations, the heliographic coordinates, the inclination angles of the loop planes, the 
loop curvature radii, the oscillation periods, and the damping times are summarized in Table 1, extracted 
from Tables I and II in Aschwanden et al. (2002), as well as one damping time from Nakariakov et al. 
(1999). 
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2.1. Parameterization of Loop Skin Depth 

Damping of oscillations and waves by resonant absorption has been studied as a mechanism for coronal 
heating. Most studies in this context are concerned with driven waves. The interest for the present paper 
is in the eigenmodes damped by resonant absorption. Hollweg and Yang (1988) derived an analytical 
expression for the decay time in planar geometry for an equilibrium model with a thin non-uniform 
boundary layer. They translated their Cartesian result to cylindrical flux tubes and were the first to point 
out that kink mode oscillations undergo fast damping. In our view Hollweg decay is a good name to refer 
to this fast damping due to resonant absorption. Goossens et al. (1992) derived analytical expressions 
for the damping rate for 1-D cylindrical flux tubes with thin non-uniform boundaries (TB) under rather 
general conditions of the equilibrium magnetic field and stationary equilibrium velocity field. Of particular 
relevance for the present study is their Eq. (77) for a static loop with a straight field. It was derived 
under the assumption that the loop is long so that the tube is thin (TT). Ruderman and Roberts (2002) 
reconsidered the problem as an initial value problem. They arrived at an analytical expression for the 
damping rate which is a special case of the corresponding equation of Goossens et al. (1992) (see their 
Eq. 56) . Ruderman and Roberts then specialized to a sinusoidal variation of density in the thin nonuniform 
layer (their Eq. 71) and obtained the corresponding decay time (their Eq. 73). Their density profile n(r) 
across a loop cross-section is parameterized by 



n(r) 



rii for r < (a — I) 



n. 



(1+X) _ (1-X) „:„ tt (2r+i-2a) 
2 2 hm 2 I 



for (a - I) < r < a (1) 



n e for r > a 



where \ = n &l n i is the density ratio of the external plasma n e to the internal plasma density rii of the loop. 
The depth I of the loop surface region might be called the "skin depth", because it characterizes the depth 
of the outer envelope layer over which the density varies. So, a is the outer loop radius, b — a — I the inner 
loop radius, and R the mean radius, which defines also the half width (wi oop /2) for the loop, 

R =^T = {^¥)= a - l 2- < 2 > 

So the density outside of the loop is n(r > a) — n e , the skin region is bound by b < r < a, and the density 
inside this skin depth is n(r < b) — rii. An example of such a density profile is shown in Fig. 2 (bottom 
panel), for an inner radius b = 0.4 Mm and an outer radius a = 3.5 Mm. 

With this parameterization, Ruderman & Roberts (2002) derive a ratio of the exponential damping 
time to to the oscillation period P (their Eq. 73), where we can replace R w a in the thin-boundary 
approximation, 

(tD_\ _ 2 /on (1 + X) „ 2 (R\ (1 + x) 

\PJmn Al) (1- X )~ (1- X )> (6) 

Ruderman and Roberts (2002) obtained the above analytical expression (and also their more general 
expression (their Equation 56)) under the assumption that the non-uniform layer is thin, meaning l/a <C 1. 
In the present paper we use the Ruderman and Roberts formula to estimate the ratio \ f° r loops with thick 
nonuniform layers. A generalization of this result for thick boundary layers thus involves two corrections, 
one for the replacement of the outer loop radius a with the mean loop radius R — a — 1/2, and a second one 
that quantifies a correction factor qxB between the thick-boundary and thin-boundary treatment using the 
mean loop radius R. Thus, the damping formula may be written in terms of R and the correction factor 
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Qtb, 

^) thic r qTB (^) t u n =qTB ^T) { ^- & 

where the correction factor qtb depends on the boundary thickness ratio (l/R) as well as on the density 
ratio x, and has been calculated numerically in Van Doorsselaere et al. (2003). For instance, for a density 
ratio of x — 1/3, the correction factor qTB varies in the range of [0.75, 1.18]. In the fully non-uniform limit 
l/R = 2, the correction value is qmil/R = 2, x = 1/3) ~ 0.75. We use the numerically calculated values 
qtb in Table 2 to predict the external plasma density. 

In a previous study we measured the oscillation periods P and damping times t D of 11 events 
(Aschwanden et al. 2002). Here we attempt to measure the loop geometry parameters a and b = a — I, and 
the density ratio x to test this theoretical model (Eq. 4) of damping by resonant absorption. The density 
contrast of the oscillating loop is just the inverse ratio x 1 = n i/ n e, which is larger than unity for every 
detectable loop. 



2.2. Loop Density Profiles 

In order to measure cross-sectional density profiles n(r) of coronal loops observed in optically-thin 
EUV or soft X-ray emission, at least four effects play a role that need to be taken into account: The 
subtraction of the background flux from the plasma in front and behind the oscillating loop along the 
line-of-sight (Fig. 2), the line-of-sight integration of the emission measure (Fig. 3), the spatial smearing 
due to the transverse motion of an oscillating loop during an exposure time (Fig. 4), and the point-spread 
function of the instrument (Fig. 5). 

We start with the background subtraction, which is simply done by inspecting cross-sectional density 
profiles, averaged over some loop segment along the loop, and by selecting the lowest flux values on both 
sides of the oscillating loop, and interpolating a linear function between both sides (Fig. 2, top panel). 
Thus, the total EUV flux per pixel across a loop cross-section is defined by, 

F(r, t) = F back (r, t) + Fi oop (r, t) , (5) 

which yields a time-averaged background flux F back (per pixel), 

F b ack =< F back (r,t) > , r left (t) < r < r right (t) (6) 
and a time-averaged loop flux Fi oop at the central axis of the loop, 

Fi 00 p = max[Fi oop (r,t)} , n eft (t) < r < r nght (t) (7) 

where the loop boundaries [ri e ft{t),r r i g ht{t)] vary as function of time t depending on the motion of the 
oscillating loop (see Fig. 6, left panel). The loop flux is related to the electron density n(z) along the 
line-of-sight z by the emission measure EM at pixel position (x, y), 

dEM(x, y, T) 



dT 



j n 2 (x,y,z,T) dz . (8) 



The observed flux Fi oop (x,y) in a given filter specified by a temperature-dependent instrumental response 
function R(T) is obtained by integrating the emission measure EM{T) with the response function R(T) 
over all temperatures T, 

F loop (x iy) =J dEM ^ R(T)dT. (9) 



- 5 - 



For an isothermal loop segment that is ncar-perpcndicular to the line-of-sight, we can obtain the radial flux 
profile F(r) by integrating the density profile n(r) specified in Eq. 1 along the line-of-sight z (Fig. 2 bottom 
and Fig.3), 



Fi oop {r) = R(T) / n 2 (r' = Vr 2 + z 2 ) dz . (10) 



The circular cross-section essentially causes a convolution of the density profile n(r) with a column depth 
Az that has a circular dependence Az oc \J (a 2 — r 2 ), so that the flux profile F(r) looks more gaussian-likc 
(Fig. 2, middle panel) than the trapezoidal shape of the density function (Fig. 2, third panel). We fit the 
(normalized) theoretical loop profile (Eq. 1) to the observed flux profiles Fi oop (r) by optimization of the 
parameters a and I, using a Powell minimization method (Press ct al. 1986). 

The measured total flux at the center (r = 0) of the loop is 

Ftotai = [njwioop + n 2 e (z back - wi oop )] R{T) , (11) 

with n, being the internal loop density wi oop the mean loop width (Eq. 2), n e being the external or average 
background density extended over a column depth zt, ac k- The mean background flux measured at the sides 
of the loop is, 

F b ack = [n 2 e z back ] R(T) , (12) 

yielding a background-subtracted flux of 

Fsubtr = F tota i - F back = [(n 2 - n 2 e )wi oop ] R{T) , (13) 

yields then the difference of the squared densities. Thus we cannot determine the internal loop density rn 
directly, but only as function of the external density n e , 



AF sub tr 

wi oop R(T) 



(14) 



In the limit of vanishing (background-subtracted) loop flux, F subtr *— ► 0, the densities inside (n,) and outside 
of the loop (n e ) become identical. The density profile fitting is performed for every time step (typically 
10-30 images) of the 11 oscillation events (e.g., see Fig. 6). The results of the best-fit parameters a, I and 
\Jn\~ n 2 , with the mean and standard deviation from averaging over all (~20-30) time steps, are given in 
Table 2. 

The data analysis procedure is illustrated in Figs. 1-8. Fig. 2 shows a fit of a cross-section profile 
n(r,ti) to the observed flux Fi oop (r,t = ti) at a single time step ti. Fig. 6 shows the fits as function of 
time ti, i = 1, n, and Fig. 7 shows the variation of the measured parameters as function of time, Fi oop (t) 
and Fback(t) (Fig. 7 top), the oscillation amplitude A(t) (Fig. 7 middle), and the cross-section parameters 
a(t), l(t), and a(t) — l(t) (Fig. 7 bottom), with the means and standard deviations indicated. Fig. 8 shows a 
summary plot of the average cross-section fits, for each of the 11 events. 



2.3. Oscillatory Motion Smearing 

We have to be aware that every TRACE image has been recorded with a finite exposure time of typically 
At exp w 5 — 10 s. Since we are measuring the loop profiles in perpendicular direction to the loop axis, the 
oscillatory motion of the transverse kink mode introduces a smearing that transforms a rectangular density 
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profile into a trapezoidal profile (Fig. 4), if not corrected. For typical oscillation speeds of v max w 10 — 100 
km s _1 (see Table 2 in Aschwanden ct al. 2002) we expect a motion of Ar = v max At exp w 50 — 1000 km, 
which corresponds to Ar « 0.1 — 2.8 TRACE pixels with a pixel size of 0.5". However, the observed loop 
radii were found to be in the range of a w 2000 — 12, 000 km, so at least an order of magnitude larger. For 
most times the actual speed is smaller at an arbitrary phase of the sinusoidal oscillations, v(t) < v max . We 
measured the actual amount of smearing for every fit and found that it amounted indeed in all cases to a 
fraction of less than 0.05-0.1 of the loop width. Therefore we neglected this effect in the fitting procedure. 



2.4. Point Spread Function 

The instrumental point-spread function contributes to some broadening and smoothing of observed 
density profiles, and thus could affect our inversion of loop density profiles from observed flux profiles. 
The point-spread function of TRACE has been investigated in 171 A EUV image fits using a blind 
iterative deconvolution (BID) procedure (Richard Nightingale, private communication). The shape of the 
TRACE point-spread function was found to have the shape of a 4-sidcd pyramid with a square-shaped 
base rotated by 45° with respect to the CCD raster (Fig. 5). The point-spread function falls off 
from a central pixel with value 1.0 to 0.33-0.36 in the next-neighbor pixels, and almost to 0.0 in the 
second-next neighbor pixels, for a pixel size of 0.5". Thus, the average full-width of the pyramid shape is 
FWHM = 2 x 0.5" (1 -0.345)/0.5 = 1.3". 

Independently, the TRACE point-spread function has also been characterized with a BID procedure by 
Golub ct al. (1999), who also found a slightly elongated shape at a position angle of 45°, with a FWHM of 
3 pixels in one direction and 2 pixels in the orthogonal direction, yielding a FWHM = 2.5 x 0.5" = 1.25", 
which is consistent with the former measurement. 

Since our measured loop widths (w = 2a — I rts 2. ..14 Mm, see Table 1) are in the average at least an 
order of magnitude larger than the FWHM of the point-spread function (FWHM w 1.3 x 0.725 = 0.94 
Mm, we neglect it in fitting of the density profiles to the flux profiles. 



2.5. Predicted External Plasma Density 

After we have measured the loop profile parameters a and I (Section 2.1 and Table 2), and using the 
measurements of the observed loop oscillation periods P and damping times to from the previous study 
(Aschwanden et al. 2002), we have only one free parameter left in the damping time expression (Eq. 4), 
namely the external-to-internal density ratio \ = n e/n,i of the oscillating loop. Because it is difficult to 
measure the ambient plasma density n e of an oscillating loop, we do not explicitly predict the damping 
time ratio tp/P based on uncertain densities n e with relation (12), but rather do it the other way around 
by using the theoretical relation (12) to predict the ambient density n e , which can then be compared with 
observational measurements. 

As Eq. (4) shows, the shortest damping time ratios occur for a loop in vacuum, i.e. for \ = 0- We list 
these minimum ratios 

in Table 2, which are calculated from the measured values of R — a — 1/2 and I and using the fully 
non-uniform approximation qrB(l/R ~ 2) w 1.0 (Van Doorssclacrc ct al. 2003). The resulting values of 
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(tD/P)min (Table 2; 7th column) reveal that they are all lower than the observed values (Table 2; right-most 
column), as expected for x > 0. This is a first successful test of the theoretical model, in the sense that all 
11 observed cases are able to provide a physical solution, namely a positive value for the density contrast, 
X>0. 

In a next step we express the density contrast xd explicitly as function of the other variables from 
Eq. (4), (where the subscript in xd indicates here that it is derived from the damping time to, instead of 
the standard definition in terms of density contrast, x — n e/ni, as defined in Eq. 1), 

XD = g^> X .-LL (>)(%). (M) 



(X + l) ' q TB 2 \RJ \P 

We predict now the density contrasts xd based on the measured ratios of damping times tjj to periods P 
((tD/P)obs in Table 2), using the fully non-uniform limit (qTB(l/R = 2) = 1.0). We find values in the range 
of xd ~ 0.3 — 0.8 (Table 3, 8th column). Using the numerically calculated correction factors qrB{l/R) 
computed with the LED A code (Van Doorsselaere et al. 2003) for the actual observed values of (l/R), 
indicated with xleda in Table 3, we see that the approximation qxB ~ 1-0 is a very good approximation 
for fully-nonuniform loops (1/R=2). 



2.6. Measurement of the External Plasma Density 

In a next step in our analysis we attempt to estimate the external plasma density n e around the 
oscillating loops from the measured background flux Fb ac k and loop position. The flux of the background 
is composed of the emission measure along the line-of-sight in front and behind the oscillating loop. We 
assume a stratified atmosphere for the spatial and temporal average of the background flux, with an 
exponential density scale height Ay corresponding to a mean temperature T. For the plasma detected in 
the TRACE 171 A passband, this mean temperature is T w 1.0 MK. Thus the vertical density profile of the 
detected coronal plasma is, 

n e (h) = noexp (- A) , At = 47 Mm (y^) , (17) 

We need to calculate the column depth of a hydrostatically stratified atmosphere along a line-of-sight as 
function of the distance d from Sun center. We define an equivalent column depth, z eq (d 7 T), as function of 
the distance d = r Q + h from Sun center, for a mean coronal temperature T, 

/oo 
n 2 e (h[z],T)dz = n 2 (T) z eq (d, T) . (18) 
-oo 

with rio(T) the coronal base density as defined in Eq. (6). From Eqs. (6-7), the following relation follows 
for this equivalent column depth (Table 3, 4th column), 



z (d,T) = f * exp(- 2[Vrf2+ N Z2 -^Ql) dz for d > r . 

Jzi *T 



(19) 



with the integration limits Z\ — — oo and z 2 = +oo for abovc-thc-limb locations (d > r Q ). Inside the disk 
(d < tq), we have only to change the integration limit z 2 to, 



z 2 (d) = -^Jr 2 Q -d 2 forO<d<r . (20) 
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The column depths z eq (d,T) are shown in Aschwandcn & Acton (2001; Fig. 3 therein) for a height range 
from disk center {d = 0) to one solar radius outside the limb (d = 2r Q ), for temperatures in the range of 
T = 1.0 — 4.0 MK. At disk center (d = 0), the equivalent column depth matches the emission measure scale 
height, which is the half density scale height (Xem = At/2). At the limb (d = r Q ), there is in principle a 
discontinuous change by a factor of 2, which, however, is difficult to measure because of the extremely high 
instrumental resolution required to resolve this jump. Above the limb, the column depth drops quickly with 
height. 

We can now relate the observed background flux Fb ac k(d), measured at a line-of-sight with distance d 
from disk center, to the emission measure EMt, ac k(d), using the instrumental response function R(T), 

Ftackid) = EM back (d)R(T) . (21) 

and determine the coronal base density uq (with Eqs. 9 and 10), 



na V R(T) Zeq {d,T) ■ (22) 

To estimate the ambient density around the oscillating loop in height h osc , we have to use the hydrostatic 
model of Eq. (6), 

n e (h = h osc ) = n cxp ( _^£f£ j ) (23) 



Obviously we need measurements of the center-limb distance d of the location of the oscillating loop 
segment, as well as an estimate of the altitude h osc of the oscillating loop segment. In the previous study 
(Aschwanden et al. 2002) we measured the heliographic longitude difference (Zi — Iq) and latitude difference 
(b\ — bo) of the midpoint of the loop baseline to disk center, the inclination angle i? of the loop plane, and 
the loop curvature radius R CU rv We list these parameters of our 11 analyzed loops in Table 1. With these 
geometrical parameters we can now determine the projected distance d of the location of the oscillating 
loop segment to disk center (Table 3, 3rd column), 

d = (r Q + h osc ) sin [(h - hf + (6i - b ) 2 ] , (24) 

and the height h osc of the oscillating loop segment above the solar surface (Table 3, 2nd column), 

h>osc — Rcurv COST? . (25) 

Inserting these parameters d and h osc into Eqs.(6)-(12) we obtain now an estimate of the ambient density 
at the height of the oscillating loops, n e (h = h osc ). The so evaluated density values n e are given in Table 
3 (5th column), along with the inferred internal densities nj (Eq. 14) and density ratios x — n &l n i- These 
density ratios can now be compared with the predicted density contrast \d (Eq. 5) from the observed 
damping times (Table 3). The uncertainties of the derived parameters were estimated according to the 
error propagation law (Appendix A). We find that density ratio is consistent with the model of resonant 
absorption within a factor of xleda/x ~ 1-2 — 3.5 (Table 3, right-most column, excluding the lowest and 
highest extreme value). 



2.7. Temperature Corrections 



In our analysis we used the temperature of T = 1.0 MK that corresponds to the peak of the TRACE 
171 A passband, in which all the analyzed oscillating loops were detected. This peak temperature is 
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certainly representative for the background plasma along the line-of-sight, because it represents an average 
over many coronal structures which are detected in a given passband with the highest probability near 
the peak temperature of the temperature sensitivity. So the temperature should not affect any derived 
parameter based on the background plasma, such as the external plasma density n e . 

What's about the temperature inside the oscillating loops. Since the FWHM of the temperature 
response function in 171 A is about AT m /T m = (1.2- 0.8)/1.0 = 0.4, a loop is detected with a probability 
of 67% in this temperature range. The peak response function we used is Rm(T = 1.0MK) = 1.1 x 10 -26 
DN s _1 cm 5 (see e.g. Fig. 12 in Aschwanden et al. 2000). If we approximate the response function with a 
gaussian curve, single-temperature plasma structures are detected with a probability of 24% at a sensitivity 
that is less than 50% of the peak response. So, statistically, in every 4th loop we may have overestimated 
the response function by a factor of > 2, which is equivalent to an underestimation of the true loop density 
by a factor of < \/2. Therefore the resulting density contrast \ = n e /rii could be a factor of 1.4 higher for 
every 4th loop. In the statistical average, however, this temperature correction is not sufficient to explain 
the average discrepancy between the density ratios, i.e., xleda/x ~ 1-2 — 3.5. 



2.8. Coronal Filling Factors 

Another not considered effect is the spatial filling factor, which affects both the plasma determination 
external and internal to the oscillating loops. Generally, if the filling factor is less than unity, density 
derivations from the emission measure (Eqs. 8, 11, 19) result into an underestimate of the density. If both 
the internal and external plasma is subject to the same filling factor, this effect would cancel out in the 
density contrast \ = ™e/«i and no correction is needed. However, we think that the oscillating fluxtubes, 
especially those with small diameters are more likely to be solidly filled with plasma than the wide bundles 
of fluxtubes, or the background corona. Inquiring the diameters of the oscillating loops we find large radii 
a > 10 Mm only for two cases (3a and 16a in Table 2), which show the same discrepancy between \ an d 
Xd as the other cases (Table 3), so a correction by a filling factor of loops cannot improve the consistency 
between data and model either. 

On the other side we can ask whether the filling factor of the background corona has an effect on 
our model. With our stratified coronal model we applied for the temporal and spatial average (Eq. 6), 
we inferred a density contrast of % = 0.30 ± 0.16. If the background corona is subject to a filling factor 
less than unity, the true ambient density around a loop could be lower or higher. A possible bias towards 
a higher value could result in active regions, where high-density concentrations are more likely around 
oscillating loops. If we consider such a filling factor bias and assume that the ambient density around an 
oscillating loop in an active region is actually higher, the density contrast value x = n el n i increases. The 
mismatch is in the average xleda/x ~ 1-2 — 3.5, which could be reconciled with correspondingly higher 
ambient densities around the oscillating loops. This higher density does not necessarily need to be plasma 
with a temperature of T w 0.8 — 1.2 MK as detected with TRACE 171 A , it could be plasma of higher 
temperature, say in the range of T w 1.2 — 2.0 MK, as many differential emission measure distributions 
inferred in active regions suggest (e.g. Brosius et al. 1996; Aschwanden & Acton 2001). Improvements 
in the determination of the mean external density n e therefore require the knowledge of the differential 
emission measure function, which demands multi-filter data. 



-10- 



2.9. Comparison with Phase Mixing Model 

The theory for damping due to resonant absorption for thin non-uniform layers predicts that the 
damping time t D is a function of the period P, the geometric ratio (R/l), and the external/internal density 
ratio x — n el n i (Eq. 4), without any free parameter, 

RA _ 2q TB RP(l + X) 

d — ; T\ \ • y 

* i (i-x) 

We test this scaling law in Fig. 9 (right frame) and find a mean ratio of t% A /t°£ s = 0.37 ± 0.15, so theory 
and observations agree within a factor of 3, where the discrepancy probably is due to the underestimation 
of the background density when measured with a narrow-band filter. 

Alternatively, we might test the scaling law for phase mixing (Heyvaerts & Priest 1983; Roberts 2000), 
which predicts the following dependence, 




where L is the loop length, l in h the scale of inhomogeneity (which we set equal to our skin depth here) , v is 
the coronal viscosity, v^ the Alfven speed inside the fluxtube, which amounts to v^ = \/2L/P for the kink 
mode in a low-/? plasma (Nakariakov & Ofman 2001). We calculate the predicted damping times with a 
standard value of the coronal viscosity, v = 4x 10 13 cm 2 s" 1 , and plot them in Fig. 9 (left frame). We find 
an average ratio of i£ M /t°£ s = 0.79 ± 0.19, which closely agrees with the observations. Thus the model of 
phase mixing cannot be excluded as alternative interpretation. 

A corresponding test of the scaling law tjj cx P for resonant absorption and to oc {LP) 2 / 3 has been 
performed in Ofman & Aschwanden (2002) that showed also a slight preference for phase mixing. The test 
here, however, is more constrained. There are three differences to the former study: (1) we do not make 
the assumption that the spatial scale of inhomogeneity k n h is proportional to the loop length L or loop 
width w; (2) The loop widths wi oop = a + b are measured here from the deconvolved density profiles and 
not from the FWHM of the flux profiles, and (3) we measure here the (outer) loop radius a and spatial 
scale of inhomogeneity I separately, which were set equal to each other in the former study. Nevertheless, 
we obtain similar results that both models are roughly consistent with the observations. 



3. CONCLUSIONS 

In this study we modeled the cross-sectional density profiles n e (r) of oscillating loops, specified by the 
outer radius a, skin depth I, internal density rii, and external density n e . These parameters allow us to test 
the theoretically predicted relation between the damping time to, oscillation period P, geometry (a, I), and 
density parameters (n e ,rii) for the damping mechanism of resonant absorption. Because we can measure 
all these observables we have no free parameters in the model and thus are able to perform a very strict 
consistency test between theory and observations. The alternative damping mechanism of phase mixing can 
be tested with these measured parameters also, but there is a free parameter, namely the viscosity, which 
cannot directly be constrained by observations to date. Our observational test yields the following results: 

1. The means and standard deviations of our measured parameters are (see Table 4): Outer loop radius 
a = 4.5 ± 3.5 Mm, loop skin depth Z = 3.9 =b 3.1 Mm, skin depth ratio l/a = 0.85 ± 0.08, internal loop 
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density nt = (1.4 ±0.7) x 10 9 cm 3 , external loop density n e — (0.36 ±0.18) x 10 9 cm 3 , density ratio 
X = n e /rii = 0.30 ± 0.16. These are the averages of 11 oscillating loops events. 

2. In a previous study we measured the corresponding oscillation periods, P = 317 ± 114 s, the 
damping times to = 574 ± 320 s, which yield a mean ratio of tn/P = 1.8 ± 0.8. According to the 
resonant damping model of Rudermann & Roberts (2002), originally derived for the thin-boundary 
approximation and now generalized for the thick-boundary approximation by Van Doorsselaere et 
al. (2003), the observed damping times constrain (under the assumption of damping by resonant 
absorption) a density ratio of xd = 0.53 ± 0.12, which is a factor of xleda/x ~ 1-2 — 3.5 higher 
than that measured from the background fluxes with the TRACE 171 A filter. It is likely that this 
discrepancy factor results from the neglected hotter plasma with T > 1.5 MK that is not detected with 
the 171 A filter. With this correction, the damping model or resonant absorption can be considered as 
a successful theory to explain the observed damping times. Alternatively, the model of phase mixing 
is also found to be consistent with the data. 

3. The damping model by resonant absorption provides a direct diagnostic of the density ratio 
Xd — n e /ni. The observed parameters of the loop cross-section profiles vary very little, the ratio 
a/l ~ 1.18 ± 0.11 varies only by ss 10% and can be neglected in the damping formula. The 
correction factor for the thick-boundary treatment can then be taken in the fully non-uniform limit, 
<Itb{1/R = 2) ;=y 1.0. Therefore we have a very simple relation that predicts the number of oscillations 
Nose — td/P as function of the density ratio xd = n e /rii, or vice versa (Eq. 16), 

N osc = % « !1±*£ , (28) 

P 7T 1 - XD 

n e irN osc — 1 . . 

Xd — ~ ° - : • 29 
Hi ttN osc + 1 

The relation is plotted in Fig. 10 and can be used as an efficient density diagnostic. 



This study provides new support for the interpretation of damping mechanism of coronal loop 
oscillations in terms of the resonant absorption process. A new effect we learned from this study is 
the sensitivity to the density contrast between the loop and the ambient plasma. In vacuum, the loop 
oscillations would be damped within a half oscillation period, tp/P ~ 0.5. However, the higher the ambient 
plasma density is, the less severe is damping by resonant absorption, so that undamped oscillations can 
only be supported if the density contrast is very little. Asking the question why only a small subset 
of all active region loops exhibit oscillations after a global triggering event, e.g. during a flare or a 
filament destabilization (Schrijver et al. 2002), the mechanism of resonant absorption provides a plausible 
explanation that oscillations are most favored in loops with little density constrast to the ambient plasma, 
while all other loops with a large density contrast are strongly damped within a half oscillation period. 
For future work to study the role of resonant absorption we recommend to model the differential emission 
measure distribution of the coronal plasma with multi-filter data to obtain a better estimate of the coronal 
background density. 
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Appendix A: Estimates of Parameter Uncertainties 



The variables of the damping time t D and period P (Table 1) have been determined in Aschwanden 
et al. (2002) without an estimate of the uncertainty. Based on multiple trials with different background 
subtractions we estimate the uncertainty of the damping time to be of order a tr> ~ i_o/3, e.g., compare the 
result of tu = 1200 s for event la) in Aschwanden ct al. (2002) versus tu = 870 s in Nakariakov et al. 
(1999). The error in the period measurement P can be neglected because repeated fitting with different 
background subtractions reproduced this value within a few percent, so the error is much smaller than the 
error of the damping time tD- Also the errors in the parameters l\ —lo,b\ — bo, Rcurv,^ and the derived 
quantities h osc (Eq. 25) and d (Eq. 24) are accurate to a few percent and thus the uncertainties can be 
neglected. The equivalent column depth z eq (Eq.19) is a theoretical quantity that has no measurement 



For all parameters directly measured in this study, Fi oop , Fback, a, and I (Table 2), we determined the 
uncertainties <Jf 1oop , &F back , c a , and 07 from the standard deviations that resulted by averaging the fits of 
all times per event, ti, i = 1, n, with typically n w 20 — 30 time steps per event. 

The uncertainties of the derived parameters rij (Eq.14) qn = to/P (Eq.4), X (Eq.6), \d (Eq.6), n 
(Eq. 22), n e (Eq. 23), and \ — n e/ni (after Eq. 1) were calculated with the error propagation law, 



error. 




(Al) 



Oq D = QD\/ (Va/a) 2 + (<7l/l) 2 , 



(A2) 
(A3) 



X^/(a tD /t D y + (a a /aY + (a l /l) 2 




(AA) 




(A5) 




(A6) 



(A7) 
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Tabic 1. Times, locations, loop geometries, oscillation periods, and damping times of 11 oscillation events 

analyzed in Aschwanden et al. (2002). 



No. 


Date and Time 


Hcliographic 


Hcliographic 


Loop 


Loop 


Oscillation 


Damping 




of Observation 


Longitude 


Latitude 


Inclination 


Curvature 


Period 


Time 






lo — Iq 


6o - 6 Q 




Rcurv 


P 


to 






(deg) 


(deg) 


(deg) 


(Mm) 


(<0 


(s) 


la) 


1998-Jul-14 1259:57 


-15.6 


-27.6 


7.0 


47.0 


261 


870 a 


lb) 


1998-Jul-14 1257:38 


-15.5 


-26.0 


19.0 


24.0 


265 


300 


Id) 


1998-Jul-14 1257:36 


-19.5 


-24.5 


-35.0 


55.0 


316 


500 


If) 


1998-Jul-14 1256:32 


-19.6 


-24.5 


-44.0 


57.0 


277 


400 


lg) 


1998-Jul-14 1302:26 


-19.2 


-22.7 


47.0 


45.0 


272 


849 


3a) 


1998-Nov-23 0635:57 


82.3 


-27.7 


-12.0 


99.0 


522 


1200 


4a) 


1999-Jul-04 0833:17 


26.0 


-27.3 


-14.0 


74.0 


435 


600 


5c) 


1999-Oct-25 0628:56 


-22.9 


-21.3 


2.0 


53.0 


143 


200 


10a) 


2001-Mar-21 0232:44 


72.6 


-3.8 


20.0 


77.0 


423 


800 


16a) 


2001-May-15 0257:00 


22.7 


-18.3 


39.0 


68.0 


185 


200 


17a) 


2001-Jun-15 0632:29 


-48.7 


-28.0 


41.0 


33.0 


396 


400 


This value of t D = 870 s was 


measured in 


Nakariakov ct al. 


(1999) and 


is also used 


in the study of Ofman & 



Aschwanden (2002). An alternative value of tu = 1200 s was determined in Aschwanden et al. (2002). The difference 
reflects a typical uncertainty in the determination of the damping time to ■ 



Table 2. Best-fit parameters of loop cross-section fits to the same 11 events specified in Table 1. 



No. 


Loop flux 
Fsubtr [DN/s] 


Background 
Fback [DN/s] 


Loop radius 

a [Mm] 


Skin depth 
I [Mm] 


Loop density 
^/n\ - n\ [10 8 cm- 3 ] 


Minimum ratio 

(^D 1 P)min 


Observed ratio 

(t D /P)obs 


la) 


9.9± 1.6 


23.2± 1.8 


3.5± 0.4 


3.3± 0.5 


15.6± 1.3 


0.35± 0.53 


3.33± 1.11 


lb) 


20.9± 5.0 


28.0± 8.1 


3.4± 0.5 


2.9± 0.5 


22.4± 2.7 


0.41± 0.26 


1.13± 0.38 


Id) 


4.0± 3.1 


18.6± 7.5 


2.6± 0.9 


2.1± 1.1 


10.8± 4.3 


0.48± 0.59 


1.58± 0.53 


If) 


0.5± 0.3 


7.3± 3.0 


2.0± 0.5 


1.6± 0.6 


4.5± 1.6 


0.49± 0.42 


1.44± 0.48 


lg) 


7.0± 2.8 


43.2± 3.7 


3.4± 0.7 


2.9± 1.0 


12.7± 2.8 


0.44± 0.55 


3.12± 1.04 


3a) 


12.7± 5.3 


34.7± 1.3 


12.4± 4.1 


10.8± 4.3 


9.1± 2.0 


0.41± 0.40 


2.30± 0.77 


4a) 


12.6± 2.2 


54.4± 12.6 


2.8± 0.5 


2.7± 0.5 


19.8± 1.7 


0.34± 0.38 


1.38± 0.46 


5c) 


13.9± 2.9 


44.8± 6.3 


2.5± 0.2 


2.3± 0.2 


21. 9± 2.3 


0.37± 0.26 


1.40± 0.47 


10a) 


25.9± 8.4 


40.5± 6.5 


4.6± 0.5 


3.6± 0.8 


20.7± 3.4 


0.48± 0.36 


1.89± 0.63 


16a) 


2.5± 0.5 


6.3± 1.5 


10.2± 2.9 


9.0± 3.3 


4.4± 0.4 


0.40± 0.59 


1.08± 0.36 


17a) 


1.2± 0.7 


53.6± 2.0 


1.6± 0.6 


1.2± 0.7 


7.1± 2.9 


0.57± 0.40 


1.01± 0.34 
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Table 3. Model parameters of external background plasma density and density contrast. 



No. 


Loop apex 


Distance 


Equivalent 


External 


Internal 


Density 


Predicted 


Numerical 


Density 


No. 


altitude 


from Center 


Column depth 


density 


density 


ratio 


ratio 


ratio 


ratio 




h oac [Mm] 


dosc/rQ 


Zeq/r e 


n e [10 8 cm" 3 ] 


rii [10 8 cm~ 3 ] 


X 


XD 


XLEDA 


Xleda/x 


la) 


46.6 


0.561 


0.040 


3.2± 0.3 


15.9± 1.3 


0.20± 0.02 


0.81± 0.27 


0.70 


3.5 


lb) 


22.7 


0.521 


0.039 


6.0± 0.7 


23.2± 2.7 


0.26± 0.04 


0.46± 0.28 


0.43 


1.6 


Id) 


45.1 


0.553 


0.040 


3.0± 1.1 


11. 2± 4.3 


0.27± 0.14 


0.53± 0.46 


0.52 


2.0 


If) 


41.0 


0.551 


0.040 


2.1± 0.7 


4.9± 1.6 


0.42± 0.19 


0.50± 0.35 


0.49 


1.2 


lg) 


30.7 


0.518 


0.039 


6.3± 1.3 


14.2± 2.8 


0.44± 0.13 


0.75± 0.28 


0.70 


1.6 


3a) 


96.8 


1.137 


0.008 


3.1± 0.6 


9.6± 2.0 


0.32± 0.09 


0.70± 0.26 


0.64 


2.0 


4a) 


71.8 


0.675 


0.044 


2.7± 0.2 


20.0± 1.7 


0.14± 0.02 


0.60± 0.37 


0.47 


3.4 


5c) 


53.0 


0.559 


0.040 


3.9± 0.4 


22.3± 2.3 


0.18± 0.03 


0.58± 0.25 


0.50 


2.8 


10a) 


72.4 


1.054 


0.109 


1.5± 0.2 


20.8± 3.4 


0.07± 0.02 


0.60± 0.27 


0.58 


8.1 


16a) 


52.8 


0.524 


0.039 


1.5± 0.1 


4.7± 0.4 


0.32± 0.04 


0.46± 0.60 


0.41 


1.3 


17a) 


24.9 


0.860 


0.061 


6.3± 1.9 


9.5± 2.9 


0.66± 0.28 


0.28± 0.36 


0.36 


0.5 
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Table 4. Means and standard deviations of measured parameters in 11 oscillation events. 



Parameter 



Mean and standard deviation 



Loop curvature radii R cur v = 
Oscillation period P = 
Damping time to = 

Observed number of oscillations tp/P = 

Predicted minimum of ratio (tjy / P) m in = 

Outer loop radius a = 

Inner loop radius a — I = 

Mean loop width wi oop = 2a — I = 

Loop skin depth I = 

Relative loop skin depth l/R = 

Loop density n; = 

External plasma density n e (T = 1MK) = 
Predicted external plasma density n e = riiXLEDA 
Density ratio \ = n e (T = 1M K)/rii = 
Predicted density ratio Xleda = le/ii = 
Prediction ratio n e /n e (T = 1MK) = Xleda/x = 



57 ± 21 Mm 
317 ±114 s 
574 ± 320 s 
1.8 ±0.8 
0.32 ± 0.05 
4.5 ±3.5 Mm 
0.6 ±0.5 Mm 
5.1 ± 3.9 Mm 
3.9 ±3.1 Mm 
1.5 ±0.2 
1.4 ±0.7 10 9 cm- 3 
0.36 ±0.18 10 9 cm" 3 
0.76 ±0.36 10 9 cm" 3 
0.30 ±0.16 
0.53 ± 0.12 
2.5 ±2.1 



- 18 - 




Fig. 1. — Oscillation event No. la on 1998-Jul-14, 12:45 UT, analyzed in Aschwanden et al. (2002) and 
in Nakariakov et al. (1999). A difference image is shown, where the transverse oscillation amplitude is 
analyzed in an area marked with a white box that is oriented perpendicular to the loop (top left panel). 
The 3D geometry of the loop is approximated with a coplanar circle (thin line in bottom panel left). The 
oscillation amplitude is decomposed into a non-oscillatory trend (top right panel) and into an oscillatory 
damped function (bottom right panel). For further details see Aschwanden et al. (2002). 
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Fig. 2. — The radial flux profile F(r) for event la (see Fig.l) is shown (top panel), as function of the cross- 
sectional radius r perpendicular to the loop and averaged along the loop within the white box shown in 
Fig.l (top left panel). A linear background to the oscillating loop is evaluated (dashed line in top panel) 
and subtracted (middle panel). A trapezoidal density function with sinusoidal boundaries (Eq.l) with outer 
radius a and inner radius a — I is shown (third panel) and fitted to the background-subtracted flux cross- 
section (middle panel) with proper line-of-sight integration across the 2D density distribution (bottom). 
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Fig. 3. — The geometry of line-of-sight (LOS) integration is shown: The 2D density distribution of a loop 
cross-section is contoured in grey-scale in the bottom left part, cross-sections along the line-of-sight are 
projected to the right side, and the radial cross-section in transverse direction (in the plane of the sky) is 
projected to the top side. 



.t=t) 



,t=y 



t=t») 



t=t 7 ) 



-> 



-> 



Fig. 4. — The effect of time smearing is demonstrated for a near-rectangular cross-section profile F(t, t = ti), 
i = 1, 7, which in superposition adds up to a more trapezoidal cross-section profile < F(r) >. 
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Fig. 5. — The point-spread function of the TRACE 171 A filter measured by a Blind Iterative Dcconvolution 
(BID) algorithm. The mesh corresponds to CCD pixels with a pixel size of 0.5" (courtesy of Richard 
Nightingale) . 
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Fig. 6. — Cross-sectional fluxes F(r, ti) for event la) are shown for time steps i=l,...,24 (left panel), averaged 
within the white box area indicated in Fig.l (top left). The flux profiles are incrementally shifted by 5 flux 
units for clarity. The background fluxes are indicated with dotted lines, the loop center location ro with a 
vertical solid line. The best-fit flux profiles F(r — ro,ti) are shown in the right panel (thick lines), based on 
the density profile n e (r) (Eq. 1) obtained by forward-fitting of the line-of-sight integration and instrumental 
response function. Note that the cross-sectional flux profile is almost invariant during this entire time interval 
of the observed observation. 
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1300:00 1310:00 1320:00 1330:00 
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-/r \ / 1 



1300:00 1310:00 1320:00 1330:00 

Time [UT] 



Fig. 7. — Physical parameters of the best-fit cross-section profiles as function of time for event la): 
Total flux Ftotaiit) = F back (t) + Fi oop (t), background flux F hack (t), and loop flux Fi oop (t) (top panel); 
oscillation amplitude A(t) (middle panel); and density profile parameters a(t) and l(t) (Eq.l) (bottom panel). 
Acceptable fits (with deviations of df / f < 10% in the flux profile) are marked with diamonds. 
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i) 1998-Jul-14 1259:57 
= 3.5 Mm 
= 3.3 Mm l\\ 



1g) 1998-Jul-14 ' 
1 =3.4 Mm 
= 2.9 Mm 



" 10a)2001-Mar-21 0232:44 
) a =4.6 Mm 
I =3.6 Mm 



)) 1998-Jul-14 1257:38 
= 3.4 Mm 
= 2.9 Mm 



■ 3a) 1998-Nov-23 0635:57 
a =12,4 Mm 
=10.8 Mm 



16a) 2001-May-15 0257:00 
1.0 a =10.2 Mm 
I =9.0 Mm 



■ 4a) 1999-Jul-04 0833:17 
i =2.8 Mm ! 
= 2.7 Mm 



" 17a) 2001-Jun-15 0632:29 

) a =1.6 Mm 

, I =1.2 Mm '■}{ 



1!) 1998-Jul-14 1256:32 
a =2.0 Mm 
I =1.6 Mm 



" 5c) 1999-0d-25 0628:56 
Ofa =2.5 Mm 
= 2.3 Mm 



Fig. 8. — Best- fit density profiles n e (r) as function of radial distance from loop center r [Mm] for all 11 
events. The dates of the events and the best-fit parameters a and I are indicated in each panel. Each profile 
is averaged from 10-30 timesteps during the time interval of the observed loop oscillation. 




Observed damping time t D [s] Observed damping time t D [s] 



Fig. 9. — Left panel: Scaling of loop oscillation damping time tjj with (a • P), which does not show a close 
correlation. For damping by phase mixing a scaling of to oc (P * l) 2 ^ is expected. Right panel: Scaling of 
loop oscillation damping time to with (R/l)P, for which a linear correlation is expected in the framework 
of resonant absorption. The linear regression fit shows a slope of 0.84 ± 0.34, which is consistent with the 
expected proportionality. 
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Fig. 10. — The damping mechanism of resonant absorption provides a density diagnostic of the density ratio 
n e /rii (of external to internal density in the oscillating loop) as function of the number of oscillation periods 
(measured by the ratio of damping time to the period, N osc = to/ P. The plot shows the prediction in the 
fully non- uniform limit (l/R = 2 and qxB = 1-0). 



